A scaling theory of 3D spinodal turbulence 
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A new scaling theory for spinodal decomposition in the inertial hydrodynamic regime is presented. 
The scaling involves three relevant length scales, the domain size, the Taylor microscale and the 
Kolmogorov dissipation scale. This allows for the presence of an inertial "energy cascade", familiar 
from theories of turbulence, and improves on earlier scaling treatments based on a single length: 
these, it is shown, cannot be reconciled with energy conservation. The new theory reconciles the i^' ^ 
scaling of the domain size, predicted by simple scaling, with the physical expectation of a saturating 
Reynolds number at late times. 

PACS numbers: 64.75+g 



A binary fluid mixture will undergo phase separation, if 
the two fluids are mutuaUy repulsive, below some critical 
temperature. Presented here is a theoretical study of the 
dynamics of spinodal decomposition in three dimensions, 
in a 50/50 mixture where the two fluids are incompress- 
ible and have the same shear viscosity, ry, and density, p. 
Starting from a completely mixed state quenched to far 
below the critical temperature, the initial separation is 
dominated by diffusion until well-defined interfaces form 
between two interlocking domains of single fluid regions. 
The ensuing late-stage coarsening is then driven by the 
interfacial tension, resisted by viscosity as the bulk fluid 
flows so as to flatten the interfaces and enlarge the do- 
main size. The coarsening can be followed through the 
average domain size, L, most commonly measured by the 
inverse first moment of the spherically averaged structure 
factor, S{k) of the difference in the concentrations of the 
two fluids, L — 2tt JS{k)dk/ JkS{k)dk. Using a simple 
(single length scale) scaling theory, Siggia [|l[ predicted 
that the domain size first shows linear growth, L ~ i, 
in the viscous hydrodynamic regime. Furukawa Q ex- 
tended this, predicting a crossover to a slower growth 
rate of t^ " in the inertial hydrodynamic regime. These 
growth rates have been observed in numerical simulation, 
see, e.g., P-Q. Linear growth has been observed experi- 
mentally P|. 

In what follows, I show that the simple scaling the- 
ory is inconsistent with energy conservation in the iner- 
tial regime; a minimal alternative is presented, based on 
three relevant length scales, which allows both force bal- 
ance and conservation of energy to be maintained at late 
times. The new scaling theory recovers the L '-^ t^'^ scal- 
ing, but gives altered scalings for other quantities (such 
as velocity gradients). This allows the physical require- 
ment, recently emphasized by Grant and Elder [^1, that 
the fluid Reynolds number should not diverge in the long 
time limit, to be reconciled with the t^/^ scaling. 

The scaling approach developed here ignores the 
anomalous scaling corrections of modern turbulence the- 
ory p9. The analysis of spinodal decomposition pre- 
sented below is thus an approximate one, but as legiti- 



mate as the Kolmogorov theory of homogeneous turbu- 
lence. As such, it is likely to be a useful tool for interpret- 
ing simulations and experiments on the inertial hydrody- 
namic regime of spinodal decomposition, and may be the 
best description available until the (simpler) problem of 
homogeneous turbulence is fully solved. 

The system can be described by the isothermal incom- 
pressible Navier-Stokes equation (NSE), 



(9v 
P-^+p(v- V)v 



T^V^v- V-P, 



(1) 



where v is the fluid velocity. (This can be represented 
by a single variable regardless of fluid composition since 
the two fluids have identical properties.) Included in the 
pressure tensor, V, is the interfacial stress which comes 
from capillary forces; the excess Laplace pressure in the 
curved interface drives the bulk fluid away from regions 
of tight curvature which eventually collapse into narrow 
"necks" that break, leading to further enlargment of the 
remaining bulk domains. Using simple scaling arguments 
p],p|pl]], the interfacial force density, V • P, can be ap- 
proximated by cr/L^, where a is the interfacial tension 
between the two fluids. The interface is assumed to re- 
main locally smooth and completely percolated through- 
out the phase separation process, with constant a. This 
is reasonable provided diffusion is rapid on the scale, ^, 
of the interfacial width (to maintain local equilibrium on 
the time scale of the interfacial motion) while being negli- 
gible over scales of the order of the domain size so it does 
not contribute to the domain growth rate. This requires 
care in simulation work |Q, but in real fluids, L ^ ^ can 
easily be arranged. 

The linear scaling is readily obtained from Eq. dl^) by 
neglecting the inertial terms on the l.h.s. and equating 
the viscous term to the interfacial one, using v ^ L and 
V ^ 1/L. When the inertial terms are no longer negligi- 
ble, reversing this argument and equating the interfacial 
force to the inertial terms, produces L ^ t'^^^. This as- 
sumes that there is only one relevant macroscopic length 
scale in the system. (Clearly there is also a microscopic 
length scale in the width of the interface, but it will be as- 



sumed, as usual, that this does not affect the macrosopic 
growth |l|.) 

A Reynolds number can be defined g| by Rcl = 
pLL/i]. Clearly for L ~ t^/^, Rcl ^ i^^^, and thus Rcl is 
predicted to grow without bound in the inertial regime. 
Recently, Grant and Elder [g| suggested that this un- 
bounded growth of the Reynolds number is unphysical, 
and there should be a further crossover of the domain 
growth rate to t^/^ (or slower) in order that Rcl remains 
finite. However, the first simulation studies to reliably 
reach the inertial regime [fj did not find evidence of any 
final growth rate slower than t^^^ . Below it is argued that 
Rbl does not quantify the inertial effects which occur in 
the bulk fluid, which are instead properly characterised 
by the Reynolds number, Rcv = p|v- Vv|/(77|V^v|), that 
is, the ratio of the nonlinear to the viscous terms in the 
NSE, and that while Re^ diverges, RCv does not. It is 
easily seen that Rcv must remain finite. Infinite Rcv im- 
plies an infinite energy density, while the fluid mixture 
starts with a finite energy density from the free energy 
difference between the mixed and separated states. 

I start from the energy equation for the system, which 
(ignoring heat flow) may be written. 



d_ 
dt 



(pw V2) + V • V(pwV2) = 

- ry (Vv)^ + V V ■ T' + {t]/ p)V^{pv^ /2), (2) 



where (Vv) = (Vv) : (Vv). Since we are not concerned 
with the local convective or diffusive energy flows, it is 
more convenient to average over the whole system and 
write the global energy balance per unit volume as. 
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(3) 



where Cin is the rate of energy transfer to the fluid 
from the interface, which can be approximated by aLjL? 
(force • velocity). 

I now show that, in the inertial hydrodynamic regime, 
the simple scaling theory is inconsistent with energy con- 
servation, by considering the behavior of the global en- 
ergy balance under the scaling predicted for this regime 
[g|. Applying the simple scaling to each term in Eq. (||), 
and replacing L by l?!"^ gives, 
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The dissipation term {—rjt~'^) clearly becomes negligible 
compared to the other two terms, as originally assumed. 
However, the kinetic energy in the fluid and the energy 
stored in the interface are both decreasing over time; thus 
energy conservation cannot possibly be maintained ex- 
cept by including the dissipation term in the energy bal- 
ance equation. Likewise, the viscous term should never 
be neglected in the NSE because it involves the highest 
order in derivatives and is thus a singular perturbation: 
however small rj is, the asymptotic physics is radically 
altered from that with t] = 0. 



I now allow for more general scaling behavior in the 
NSE and global energy balance equations by introduc- 
ing two new lengths, iv smd L^ji, with associated scal- 
ing exponents a' and a" , for the velocity first and sec- 
ond derivatives respectively, and allowing the velocity to 
scale as f^. Terms associated with the interface are still 
assumed to have scaling determined by the domain size. 
The scaling quantities are defined as follows, 

domain size: L ^ t" 

fluid velocity: v ^ t^ 

velocity first derivative: Vv ^ v/Ly "^ t^~" 

velocity second derivative: V^v ~ v/Ly2 ^ t^^'^" . 

This supposes that the interface remains smooth, but 
allows nontrivial structure in the fluid velocity at smaller 
scales. Using these scalings to write the NSE, Eq. (Q), 
and energy balance equation, Eq. (g), as powers of t 
gives, respectively, 

NSE: pl3t^-^ + pt^^-"' - Tyi'^-^"" + at-^" (5) 

energy: p/3i2/5-i ^ -ryt^'^-^o' + ^^-a-i^ (g) 

where the prefactors have been left in to facilitate iden- 
tification of the terms. This more general scaling ignores 
the anomalous scaling corrections of modern turbulence 
theory |TQ] , but is sufficient to make progress on the prob- 
lem of spinodal decomposition. The local energy equa- 
tion, Eq. (0), becomes. 
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Solutions for a, a', a" and (3 in Eqs. (|^) and (||) will 
later be checked in Eq. (R) to confirm that there are 
no discrepancies predicted for the behavior of the local 
energy flows within this three-length-scaling analysis. 

As already pointed out, dissipation must remain sig- 
nificant, so we look for a three-way balance between the 
terms in the energy balance equation, Eq. (@) , giving for 
the exponents. 



2P-l = 2f3-2a' = -a- 1. 



(8) 



This gives a' = 1/2, and f3 = —a/2. Substituting these 
back into the NSE, Eq. (||), gives. 
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^^-a/2-2«" +crt-2". 



(9) 



There is no solution to this with all four terms having the 
same exponent; solutions can instead be found by balanc- 
ing the terms off in pairs. Numbering the terms 1-4 from 
left to right, the pairing that gives the inertial regime 
scaling is term 1 with term 4, and term 2 with term 3, 
giving a = 2/3, and a" = 5/12, with /? = -a/2 = -1/3. 
The four terms in the NSE thus scale as follows. 
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The physical interpretation of this is that the moving 
interface gives rise to large scale velocity motion via 
pdv/dt. The nonlinear term, pv • Vv, then transfers 
the energy from large length scales to small length scales 
where it is finally removed by dissipation. This is the fa- 
miliar "energy cascade" of turbulence theory. Note that 
scaling arguments only predict that the paired terms bal- 
ance approximately, so the nonlinear term, pv ■ Vv, is, in 
general, larger than the viscous term, r;V^v, due to the 
energy "in transit" from large to small scales, contained 
in a series of eddies of decreasing size. The transverse 
components of pdv/dt are thus also larger than required 
to balance the interfacial force, to account for the ro- 
tational motion from pv ■ Vv. The length scales asso- 
ciated with V and V^ both grow more slowly than L, 
with iv ~ t^^^, and Ly^ ~ t^^^^, so there is an increas- 
ing separation of length scales within the system. The 
dissipation is thus decoupled from the interfacial energy 
input, and no longer affects the domain growth rate. The 
Reynolds number, defined as the ratio of the nonlinear to 
viscous terms. Rev = p\'v ■ Vvl/(77|V^v|), remains finite 
(satisfying the physical demand of Grant and Elder M), 
while the domain size grows as t'^'^ (contrary to their de- 
duction that a < ^). With this scaling, the local energy 
equation, Eq. (|7|), becomes. 



P^-5/3 



pt 



-3/2 



.^t-'/3 + ^t-5/3 + !lt-V\ (11) 
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The local convective and diffusive terms are dominant, 
and balance each other, representing the energy moved 
around by the turbulent fluid flow. 

These results are summarised in Table |l alongside the 
predictions of the simple scaling theory, for comparison. 
(Results for the viscous hydrodynamic regime are also 
shown.) In particular, notice that in the new scaling 
theory for the inertial regime, the lengths Ly and L^2 
have the same scaling as A = (5?7(w^)/e)^/^, the Taylor 
microscale, and Ad = 27r(77^/jO'^e)^/'*, the Kolmogorov dis- 
sipation scale respectively. The Taylor microscale char- 
acterises the length scales in a turbulent fluid at which 
dissipation becomes significant, while the Kolmogorov 
dissipation scale marks the small-scale end of the dis- 
sipation range M]. Thus, although these new scaling 
results have been obtained without specific input from 
turbulence theory, the extra length scales for the velocity 
derivatives turn out to (within prefactors) coincide with 
key characteristic quantities in turbulence phenomenol- 
ogy, providing strong support for the new theory. 

It is easily shown that the familiar linear scaling of the 
viscous regime M is consistent with energy conservation 
by substituting a — a' = a" = 1 and /3 = into Eqs. (0) 
and d), giving. 
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respectively. In the energy equation, Eq. (13), there is a 
simple balance between energy input and dissipation, but 



in the NSE, Eq. (p2|), the nonlinear term on the l.h.s., 
which was assumed to be negligible in the original simple 
scaling argument, appears to be decaying more slowly 
[t^^] than the r.h.s. terms (t~^). Recalling that the 
viscous regime is not the long time asymptotic regime, 
and only expected to hold for times earlier than some 
crossover time (before which the nonlinear term will be 
smaller than the other two terms), this apparent difficulty 
is eliminated. 

There are a few further potential solutions to the ex- 
ponents in Eqs. (0) and (O), which I will now discuss 
briefly. All can be eliminated on physical grounds. In 
the inertial regime, where the NSE terms are balanced in 
pairs, the other possible pairings of terms must checked. 
It is not possible to balance term 1 with term 2 and term 
3 with term 4 because the nonlinear term (term 2) is a 
force perpendicular to v (this is obvious when it is writ- 
ten in the alternative form, — pv x V x v) and therefore 
it cannot change the magnitude of v, impying (3 = Q. 
However, term 1 is proportional to (3, so is only non-zero 
if /3 ^ 0. Matching term 1 with term 3, and term 2 with 
term 4 gives a solution for the exponents of, a = 1/2, 
a' = 1/2, a" = 1/2 and (3 = -1/4. Physically, this so- 
lution has just one length scale in the system, but the 
velocity is decoupled from the interface, v ~ t~^^\ while 
L{t) ~ t^^/^. In the local energy equation, Eq. (0), this 
solution becomes, 

_ ^i-3/2 + i-5/4 ^ _ ^-3/2 ^ ^^-3/2 ^ ^^-3/2^ (^4) 

The convective term, pt^^'"^, will eventually come to 
dominate over all the other terms, thus it cannot rep- 
resent an asymptotic solution for late times. There is 
one further solution for /3 = 0, obtained by solving Eqs. 
(H) and (P) for the remaining exponents, a, a' , and a" , 
giving, a = 1/3, a' — 2/3 and a" = 1/3. The length 
scale for velocity gradients is related to the domain size 
by iv ^ i^(0; suggesting that the nonlinear term is 
mixing on scales larger than the domain size. Physically, 
this could correspond to a "turbulent remixing" regime 
as suggested by Grant and Elder p. However, the re- 
sulting breakup of the interface is liable to invalidate the 
assumptions made in deriving the scaling approximation 
for the interfacial force. Thus, although this appears to 
be a consistent solution to the NSE and energy balance 
equations, it seems unlikely that the system could, in 
fact, ever achieve such a scaling. 

In summary, I have obtained consistent scaling behav- 
ior for domain growth in the spinodal decomposition of 
a symmetric binary fluid mixture by including two extra 
macroscopic length scales for the velocity derivatives to 
enable energy balance to be satisfied. The velocity itself 
is stfll found to scale as L{T). In a system restricted 
to derivatives up to second order (the NSE is, itself, an 
approximation based on separation of macroscopic and 
molecular length and time scales), one length scale per 
derivative would seem to be a reasonable maximum for 
the purpose of simple scaling arguments, and it does not 



seem to be possible to satisfy the NSE and energy equa- 
tions using fewer. 

With this new scahng, the inertial regime scahng of 
L ~ t^/^ is maintained while the fluid Reynolds num- 
ber, Rev, measured as the actual ratio of the nonlinear 
to the viscous terms in the NSE, remains finite. The 
Reynolds number obtained from the domain size, Re^, is 
not a good estimate of Rev in the inertial regime so it 
is not physically significant that Re^ continues to grow 
without bound. The key point for the scaling behavior 
of L is that the driving force from the interface balances 
against the acceleration term alone, so is decoupled from 
what happens in the remainder of the fluid motion. The 
nonlinear and viscous terms are free to find their own 
balance independent of the driving force, provided that 
they do, eventually, remove energy from the system. 

There is, as yet, no simulation or experimental work 
that tests this new scaling theory. The simulation work 
reported in [0,|4| , while showing hints of possible differ- 
ent scaling for the NSE terms, was not able to probe far 
enough into the inertial regime to provide a significant 
test. Using Lq = V^ /{p<^) and tg = if /{pa'^) to provide 
non-dimensional length (L/Lq) and time (t/to), Ref. 
suggests that L/Lq > 10^, t/to > 10^ will be required, 
at least an order of magnitude beyond present simula- 
tion capabilities. Experimental work, though more feasi- 
ble, will not be easy to perform. Using values for water, 
rj = 10^'^ Kg m^^ s^^, and for water-paraffin, a = 2.4 
10^'^ Kg s^^, with the projection from simulation work 
of L/Lq > 10^ to estimate the interfacial force (a IL?\ 
the density matching of the two mutually repulsive fiuids 
must be better than one part in 10^ for terrestrial grav- 
itational effects to remain small in the required regime 
(i ~ 4 cm) . The micro-gravity environment of the space 
shuttle would ease this constraint p5[ |. 
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VISCOUS 

regime 



inertial regime 



simple scaling 



new scaling 



domain size L 

length for fluid velocity first derivative Lv 

length for fluid velocity second derivative L-^2 

fluid velocity v 
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1/2 
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-1/3 
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-4/3 


■4/3 
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1/3 


1/3 





_2 


-5/3 


2/3 


1/2 


1/2 


5/12 



NSE 
terms 



pd\ /dt 
pv • Vv 



= 

0/-1 
-2 
-2 



Reynolds number from interface 
Reynolds number in fluid 
dissipation rate 
Taylor microscale 
Kolmogorov dissipation scale 



Rol = {p/ri)LL 

Rcv = pv • Vv/r;V^v 

e = 7)(Vv)2 

A = (57?(«2)/e)^/2 

Ad = 2^(',?Vp3e)i/4 



1 

= 

-2 

1 

1/2 



TABLE I. Summary of predicted scaling exponents for the viscous and inertial regimes. The new theory has the same 
predictions for the viscous regime as the simple theory, apart from the NSE term pv ■ Vv. Entries are powers of time, t; an 
entry of indicates the quantity is constant, while an entry of — indicates the quantity is assumed to be zero in the viscous 
approximation. Bold entries indicate new scaling predictions that differ from the simple theory. 



